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The decay mode B —> K*l + £~ results in the measurement of a large number of related observables 
by studying the angular distribution of the decay products and is regarded as a sensitive probe of 
physics beyond the standard model (SM). Recently, LHCb has measured several of these observ¬ 
ables using 3fb _1 data, as a binned function of q 2 , the dilepton invariant mass squared. We show 
how data can be used without any approximations to extract theoretical parameters describing the 
decay and to obtain a relation amongst observables within the SM. We find three kinds of signifi¬ 
cant disagreement between theoretical expectations and values obtained by fits. The values of the 
form factors obtained from experimental data show significant discrepancies when compared with 
theoretical expectations in several q 2 bins. We emphasize that this discrepancy cannot arise com¬ 
pletely due to resonances and non-factorizable contributions from charm loops. Further, a relation 
between form factors expected to hold at large q 2 is very significantly violated. Finally, the relation 
between observables also indicates some deviations in the forward-backward asymmetry in the same 
q 2 regions. These discrepancies are possible evidence of physics beyond the SM. 

PACS numbers: fl.30.Er,13.25.Hw, 12.60.-i 


I. INTRODUCTION 

The rare decay B —$■ K*£ + £~ involves a b —> s fla¬ 
vor changing loop induced transition at the quark level 
making it attractive mode to search for physics be¬ 
yond the standard model (SM). Indirect searches for new 
physics (NP) involving loop processes require a compar¬ 
ison of theoretical estimates with experimental obser¬ 
vations. The theoretical estimates thus need to be ex¬ 
tremely reliable in order to make a conclusive claim on 
the existence or non-existence of NP. Fortunately, sig¬ 
nificant progress has been made in understanding the 
hadronic effects involved in the decay B —i K*£ + £~. The 
mode B —> K*£ + £~ is also of special interest as it allows 
for the measurement of several observables using the an¬ 
gular distribution [1]. The large number of observables 
depend on theoretical parameters that describe this de¬ 
cay. In this paper we show how some of the parameters 
can be extracted directly from LHCb measurements al¬ 
lowing us to verify our theoretical understanding. Any 
discrepancy observed must be attributed either to a fail¬ 
ure of our understanding of hadronic effects or to the 
existence of NP. We also test the relation between ob¬ 
servables that provides another clean test for NP. 

The differential decay distribution [1-3] of B —» 
K*£ + £~ results in the measurement of at least nine ob¬ 
servables using the angular distribution, as a function of 
q 2 the dilepton invariant mass squared. These observ¬ 
ables are commonly chosen to be the differential decay 
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rate with respect to q 2 , two independent helicity fractions 
that describe the decay, the three asymmetries that de¬ 
scribe the real part of the interference between different 
helicity amplitudes and three asymmetries that describe 
the imaginary part of the interference. 

Recently LHCb [4] has reported measurements of all 
these observables that have been averaged in eight q 2 
bins. A lot of studies on this decay mode are widely 
discussed in literature [5, 6]. We use the LHCb data to 
obtain estimates of hadronic form factors that describe 
the decay. Previously some of the form factors have been 
determined [7] using lfb -1 of LHCb data. We emphasize 
that our approach does not involve evaluating the decay 
amplitude in terms of theoretically estimated parameters. 
Instead we start with the most general parametric form 
of the amplitude without any hadronic approximations 
within the SM (see Eq. (5) below). Experimental data 
alone is used to fit the theoretical parameters introduced 
in the parametric amplitude. These experimentally fit¬ 
ted theoretical parameters are simply compared to the 
estimates by other authors [10, 14] which are widely re¬ 
garded as the state of the art. The values of form factors 
obtained from experimental data show significant dis¬ 
crepancy when compared with theoretical expectations 
in several q 2 bins. 

In addition to extracting form factors from data, this 
mode also allows a relation among observables that can 
provide a clean signal [2, 3] of NP. We find that the mea¬ 
surements do not satisfy the expected relation between 
the observables in the same q 2 domains where the fitted 
form factors also show a large discrepancy with the theo¬ 
retical estimates. The simultaneous observation of these 
discrepancies points to possible evidence of NP. 

The paper is organized as follows. In Sec. II we de- 
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scribe the theoretical framework developed to write the 
most general parametric form of the amplitude and cast 
the observables in a form where hadronic paremeters can 
be obtainable from data. The relation among observ¬ 
ables are also derived here. A numerical analysis is pre¬ 
sented in Sec. Ill which contains two subsections. The 
Sub-sec. Ill A gives elaborate description of extraction of 
form factors using LHCb measurements, whereas, the va¬ 
lidity of the relations derived assuming SM are examined 
in Sub-sec. IIIB with experimental data. In Sec. IV we 
summarize the important results obtained in this paper. 
Appendix. A and B estimate the complex part the ampli¬ 
tude and the systematic uncertainty arising mainly due 
to bin average effect of the observables, respectively. 


II. OBSERVABLES AND THEORETICAL 
FRAMEWORK 


In this section we briefly discuss the theoretical frame¬ 
work derived to take into account all possible contribu¬ 
tions within SM for the decay B —> K*£ + £~. We start 
with the observables as defined in Ref. [2] to be the F B , 
F±, A 4 , Ag, ApBi Ay, A 8 , Ag and dY/dq 2 = Yf. The 
observables F±, A 4 , Ag, A FB , Ay, A§ and Ag are related 
to the CP averaged observables S 3 , S 4 , Sg, Apg C!> , St, Sg 
and Sg measured by LHCb [4, 8] as follows: 


F\ = 


1 ~F l + 2 S 3 


A 4 — S 4 , 

n 

3 „ 2 


Ag — 


Afb = —Apg Cb , Ay = — St , A 8 =— Sg, A 9 = —Sg. (1) 


7T 


27T 


It may be remarked that LHCb collaboration observes 
a local tension with some observables based on the 
hadronic estimates of Refs. [10, 11]. 

We begin by assuming the massless lepton limit but 
generalize to include the lepton mass. The corrections 
due to the mass of the leptons are easily taken into ac¬ 
count [2]. In the massless lepton limit the decay is de¬ 
scribed in terms of six transversity amplitudes which can 
be written in the most general form as, 


A L X R = C X L,R F X ~G X = T C W )F\ - Q x . (2) 

This form of the amplitude [2] is the most general para¬ 
metric form of SM amplitude for B — > K*£ + £~ de¬ 
cay that comprehensively takes into account all con¬ 
tributions up to 0(Gf) within it. The form includes 
all short-distance and long-distance effects, factorizable 
and nonfactorizable contributions and resonance contri¬ 
butions. In Eq. (2) Cg and Cio are Wilson coefficients 
with Cg being the redefined “effective” Wilson coefficient 
defined [2, 12] such that 


C A =Cg + ACg fac) (q 2 ) + AC A ’ (non - fac) ( q 2 ) (3) 


where ACg ac \q 2 ), ACg’ < ' non ~ fac \q 2 ) correspond to fac¬ 
torizable and soft gluon non-factorizable contributions. 


Strong interaction effects coming from electromagnetic 
corrections to hadronic operators do not affect Cio- 
The form factors F\ and Q x introduced in Eq. (2) can 
be related to the conventional form factors describing the 
decay as shown in the appendix of Ref. [2]. The form- 
factors F\ are of particular interest here as we show that 
they can be extracted directly from data. The F\ can be 
related to the well known form-factors V, A 4 and A 12 by 
comparing with [10]: 


V{q 2 ) 


•Fl =Nx/2J\(rn 2 B ,m 2 K ,,q 2 )- 

V D m B +rn K * 

F\\ = - Ny/2(m B +m K *)A 1 (q 2 ), 

T -~ N 

*r n — 


Vo 2 


8 m B m K *A 12 (q ). 


(4a) 

(4b) 

(4c) 


It should be noted that .Fa’s and Cio are completely 
real in the SM, with all imaginary contributions to the 
amplitude arising only from the imaginary part of com¬ 
plex Cg and g x terms. Thus with the introduction of two 
variables r\ and e x the amplitude A^’ R in Eq. (2) can 
be rewritten as, 

A L C R = {TC w ~r x )Fx+iex, (5) 


where, 


r A = - Re(Cg), 

e A = Im(Cg)F'A - Im(0 A )- 


( 6 ) 

(7) 


The imaginary contributions arise mostly from reso¬ 
nant long-distance contributions, which can be removed 
by studying only those q 2 regions where no resonances 
can contribute. In practice this means the removal 
of charmonium resonance regions from the whole q 2 
range. LHCb 3fb _1 measurements [4] conservatively 
exclude the resonance region. Moreover, the contribu¬ 
tions from imaginary parts are bounded directly from the 
LHCb measurements and the bin average values of the 
e x s are found to be very small as shown in Appendix. A. 
Hence for now we are neglecting the e x ’s and will address 
it’s contribution in the numerical analysis. 

It is convenient to define Pi and P 2 as, 


Pi = 




( 8 ) 


The observables F±, F B , A F b, A 5 and A 4 can be writ¬ 


ten [2] as 

F±=ul + 2C (9) 

F l P 2 2 = u 2 + 2C (10) 

^FB = ^2 ( M ll ^ U J-) (11) 

A 4= (2C±n 0 n||) (13) 

TTri V 2 
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FIG. 1. (color online). The allowed region for Pi versus P 2 plane. The innermost yellow (lightest), the middle one orange 
(light) and outer most red (dark) contours represent lcr, 3cr and 5<r regions, respectively. The theoretically estimated values 
using Ref. [10] for q 2 < 8 GeV 2 and Ref. [14] for q 2 > 11 GeV 2 are shown as points with error bars. The light blue bands denote 
exact solutions for the SM observables including charmonium resonances from Ref. [17] parametrization and are shown only 
for the relevant q 2 bins. In most cases, there is reasonable agreement between the theoretical values and those obtained from 
data. However, for the ranges 0.1 < q 2 < 0.98GeV 2 , 6 < q 2 < 8GeV 2 , 11.0 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 there are 
significant disagreements. 


where, 


r= ncfo 

r f ’ 
u 2 2 Fjrj 


2 F‘{ 




(14) 

(15) 


u\ is always taken to be positive and the sign ambiguities 
introduced in Eqs. (11)-(13) ensure that we can make this 
assumption. The five observables F±, Fl, Afb, A 5 and 


A 4 have been expressed above in terms of five parameters 
Pi, P- 2 , C, and u_l. The other three observables A 7 , A s 
and Ag have already been used to solve for the three 
values which are presented in Ref. [2], It may be noted 
that since Fj = 1 — F^ — F±, un is not independent and is 
related to the other parameters by, rty = P\ (l — P^ 2 (rto + 
2C)-K+2C))-2C. 

In Refs. [2, 3] a relation depending on observables in¬ 
cluding all possible effects within SM which was derived 
as, 


y 4 (f l +F|| +V2nA 4 ) Fl - y (A FB + v^ 4 5 ) “ = ^Fx-y^ + yjAF L F ± ~'^-Al (16) 


This equation can be used to express any of the observ¬ 
ables in terms of the others. A solution for A 4 [3] is 


8A 5 A F b 

Aa = — ^ r 


4Fl] F ± -fA 2 BX /4F L F ± -%At 


9 ^5 


9ttFl ' 2 v / 2ttF l 

Whereas, the solution for A 5 and Afb are given by, 


■ (17) 


TTAt^FB 

2Fn 


3j4F||Fx-^A| m /2F||F L -7r2A2 


8 F„ 


( 18 ) 


4 nJUAs , ^4F L F ± -fA^2F l{ F L -^Al 

Afb =~fT ± --' ( 9) 


It may noted that Eqs. (17), (18) and (19) depend only 
on observables and not on any theoretical parameters and 
thus provides an exact test of the gauge structure of SM 
with experimental measurements. 
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FIG. 2. (color online). The allowed region for Pi versus C, plane. The color code is same as Fig. 1. The theoretically 
estimated values from Ref. [10, 14] are shown as points with error bars. The Pi and £ values significantly disagree for ranges 
6 < q 2 < 8GeV 2 , 11.0 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 , similar to the values of Pi and P 2 shown in Fig. 1. 


III. NEW PHYSICS ANALYSIS 

In this section we discuss the detailed numerical anal¬ 
ysis using 3fb _1 of LHCb data [4]. It contains two differ¬ 
ent parts, at first we show how the experimental data can 
be used to extract out the form factors which are involved 
in this decay mode. Secondly we present the consistency 
of data to test the relation among observables derived 
relying only on the gauge structure of SM. 


A. Form Factor Extraction 


We demonstrate the technique to extract out the 
hadronic parameters by including complex contributions 
of the amplitudes and considering systematic uncertainty 
arising mainly due to bin average effect. 

It is shown in Ref. [2] that £a’s contribute to the he- 
licity fractions F\ and asymmetry A 4 . We refer to Ap¬ 
pendix. A for thorough details of evaluation of the com¬ 
plex part of the amplitudes. Using the bin average central 
values of E\/ with ±ler errors from Table. Ill we can 
numerically separate out the complex contributions from 
experimental measured values of the observables. We cal¬ 
culate the central value with ±1 <7 error of the modified 
observables Fpp' and Af x/ given by, 


^ex/ _ JP eyi 


A 


ex/ 

4 


_ jex 
— s±4 



2v / 2£o£|| 

tt!) 


( 20 ) 

( 21 ) 


which enter in the % 2 definition Eq. (23) below. It en¬ 
ables us to take into account the complex corrections in 


our analysis and extract out the variables Pi, P 2 , (, uq 
and u _l (which only deal with the real part of amplitude) 
from experimental measurements of the observables ac¬ 
curately. 

It should be noted that Eqs. (9) - (13) are valid for 
each q 2 point. However, experiments can provide bin in¬ 
tegrated values of observables over a certain q 2 intervals. 
Thus a x 2 fit with bin average values of the observables 
may lead to a biased conclusion. To avoid this issue we 
have added systematic uncertainties for each observables 
due to bin average effect with the introduction of new 
parameter /3, where the change in each observable O is 
given by, 


O^O + PO 8 . (22) 

O s is the maximum shift for each observables with a best 
fitted q 2 function to 14 bin LHCb data [4] within the 
concerned bin interval. The precise determination of O s 
is described in Appendix. B. Therefore the % 2 definition 
is 



Ff* -F ± - PF f 


A F| x/ 


A| x/ -A 4 - pA% 


A Af' 


Ar-Al-pAf 


2A| X AA| X 


/ Fg*' - F l - p Fj 
^ A F£ x ' 

A 2 F b - a 2 fb PA 2 

9 A ex A 4 ex 


(23) 


where Apg and A| x indicate experimental central val¬ 
ues of the observables Apu and A 5 with ±ler errors as 
AA|g and AA| X , respectively. Similarly F™', Ffp' and 
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Af' are the central values of the modified observables 
defined in Eqs. (20) and (21) and AFf*', AFff', AA| x/ 
are ±1 ct uncertainties in it. The systematic uncertain¬ 
ties added for each observables are denoted as F±, F[, 
Af, A|g, A| s and these values are quoted in Table. IV 
of Appendix. B. The observables F ±, Fl, A 4 , A| b and 
A l are evaluated in terms of the five parameters Pi, P 2 , 
C, u 0 and u± using Eqs. (9) - (13). Considering the 
inverse of the covariance matrix the error ellipsoids are 
constructed for all the eight bins corresponding to the q 2 
values in the range (0.1 — 0.98) GeV 2 , (1.1 — 2.5) GeV 2 , 
(2.5 - 4.0) GeV 2 , (4 - 6) GeV 2 , (6 - 8) GeV 2 , (11.0 - 
12.5) GeV 2 , (15 - 17) GeV 2 and (17 - 19) GeV 2 . It can 
be seen that /3 is treated as a nuisance parameter with 
values 0±1. The y 2 function is minimized w.r.t six pa¬ 
rameters Pi, Po, C, uo, u_ l and /? and the contours shown 
in Fig. 1 and Fig. 2 are the allowed regions in the corre¬ 
sponding planes. The minimum values of the y 2 function 
for first to eighth bins are 6.9 x 10~ 9 , 3.4 x 10 —10 , 0.055, 
8.6 x 10“ 30 , 1.094, 0.538, 0.218 and 0.044, respectively. 
The best fitted values with Tier errors of the parame¬ 
ter /3 for all eight bins are 7.4 x 10 -5 ± 0.015, 1.6 x 
10~ 5 ±0.020, 0.153±0.011, 1.0 x 10" 17 ±0.005, 0.736± 
0.020, 0.251 ± 0.003, 0.261 ± 0.001 and 0.161 ± 0.012, 
respectively. 

The contours corresponding to la, 3 a and 5a permit¬ 
ted regions for Pi versus P 2 plane are presented in Fig. 1. 
These contours are compared with the estimated values 
of Pi and P 2 using Ref. [10] for q 2 < 8 GeV 2 and Ref. [14] 
for q 2 > 11 GeV 2 . The center black point denotes the 
best fit point by minimizing the chi-square function de¬ 
fined in Eq. (23). In most cases reasonable agreement is 
found between theoretical values of Pi and P 2 and their 
values obtained from data. However, there are some sig¬ 
nificant disagreements. The values of form factor ratio 
P 2 differ by 9a in the 0.1 <q 2 < 0.98 GeV 2 bin. It may 
be noted that this region in q 2 is highly affected by finite 
lepton mass and hence the large discrepancy may not ac¬ 
curately reflect the significance due to the unaccounted 
lepton mass correction systematics. Significant devia¬ 
tions are also found for the three bins 6 < q 2 < 8 GeV 2 , 
11.0 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 where Pi 
(P- 2 ) differ by 4.2a (0.8a), 5.2a (4.8a) and 5.5a (5.3a), 
respectively. The light blue bands denote exact solutions 
for the SM observables including charmonium resonances 
from Ref. [17] parametrization and are shown only for the 
relevant q 2 bins. The detailed analysis of resonance effect 
will be discussed later in this section. 

In Fig. 2 contours similar to Fig. 1, but corresponding 
to Pi versus £ permitted regions are presented for ler, 
3cr and 5cr confidence level regions. These contours are 
similarly compared with the estimated values of P 1 and £ 
using Refs. [10, 14] and assuming the theoretical estimate 
of C 10 [15]. Data shows consistency with theoretical val¬ 
ues of Pi and ( in most cases except for the two bins 
11.0 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 where 
C disagrees by 2.8(7 and 1.7a respectively. The best fit 
value of C with ±1 a error obtained from the fit can be 


used to calculate the form factor F± using Eq. (14). 


q 2 range in GeV 2 

V(q 2 ) 

Ai(d 2 ) 

AMi 2 ) 

0.1 < q 2 < 0.98 

0.677 ±0.092 
(3.05<r) 

0.570 ±0.077 
(3.40a) 

0.246 ± 0.034 
(0.88a) 

1.1 <q 2 < 2.5 

0.625 ±0.071 
(2.78<r) 

0.409 ± 0.046 
(2.00a) 

0.326 ±0.047 
(0.69a) 

2.5 <q 2 < 4.0 

0.230 ±0.150 
(1.36a) 

0.180 ±0.118 
(1.09a) 

0.214 ±0.149 
(0.81a) 

4.0 < q 2 < 6.0 

0.552 ±0.043 
(l.07cr) 

0.400 ± 0.032 
(l.69a) 

0.359 ±0.041 
(1.09a) 

6.0 < q 2 < 8.0 

0.485 ± 0.045 
(1.27 a) 

0.598 ±0.073 
(3.18a) 

0.252 ±0.025 
(1.78a) 

11.0 < q 2 < 12.5 

0.166 ±0.018 
(5.64a) 

0.560 ±0.065 
(1.76a) 

0.450 ± 0.054 
(1.81a) 

15.0 < q 2 < 17.0 

0.828 ±0.120 
(2.79a) 

0.649 ± 0.098 
(1.38a) 

0.496 ± 0.074 
(1.51a) 

17.0 < q 2 < 19.0 

1.813 ±0.436 
(0.78a) 

0.698 ±0.171 
(0.80a) 

0.461 ±0.112 
(0.91a) 


TABLE I. The form factor values obtained from fit to 3fb _1 
of LHCb data [4] . Round brackets indicate the standard devi¬ 
ation between fitted values and theoretical estimates [10, 14]. 
Significant discrepancies are found for V and Ai in several q 2 
region. 


Finally the form factor V can be evaluated using 
Eq. (4a) and the value of F± obtained. Since the recent 
3 fb -1 of LHCb result [4] does not provide branching frac¬ 
tion measurement for the entire q 2 region we assume the 
theoretical values of 1} [10, 14] in addition to C 10 [15]. 
The form factors 7F\\ and Tg can then be determined from 
the fits to Pi and P 2 respectively, using Eq. (8). Thus 
the conventional form factors A 4 and A 42 can easily be 
estimated with the relation given in Eqs. (4b) and (4c). 
In Table. I we list the best fit values with the ler un¬ 
certainties for the three form factors V(q 2 ), A 1 {q 2 ) and 
A 42 (i q 2 ) for all the eight q 2 intervals. We also present the 
standard deviation of the fit compared to the theoretical 
estimate from Refs. [10, 14]. While sizable discrepancy 
is seen for all the form factors especially in the regions 
q 2 < 2.5 GeV 2 and q 2 > 6 GeV 2 . It is interesting to note 
the very significant discrepancy is observed in the values 
of form factors V and Ai in bins 0.1 < q 2 < 0.98 GeV 2 , 
1.1 < q 2 < 2.5 GeV 2 , 6 < q 2 < 8 GeV 2 , 11.0 < q 2 < 
12.5 GeV 2 and 15 < q 2 < 17 GeV 2 . The lattice estimate 
of the form factors currently does not include finite K* 
width. This implies, that the significance of the devia¬ 
tions can be lower if one includes the unaccounted sys¬ 
tematics due to the finite K* width. We point out that 
previous attempts to incorporate resonance contributions 
in theory has been done by parametrically taking its ef¬ 
fect in the Wilson coefficient Cg [16, 17]. However the ac¬ 
curacy of the form of resonance parametrization does not 
alter our determination of form factors since, our analysis 
is independent of Cg estimates. Cg contributes only to 
u\ s and the ratios of form factors Pi and P 2 do not get 
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FIG. 3. (color online) Illustrative plots for bin average and resonance effects in the solutions for Pi (left panel) and P 2 (right 
panel). The SM observables are assumed from lattice form factors [14]. The black ‘stars’ denote the solutions obtained at 
seven different points in q 2 for the corresponding parameters in each plot. The black central curve with gray band is the form 
factor estimate (mean with ±1 <t error) of Pi and P 2 . The blue error bars are the solutions for Pi and P 2 using the bin average 
values of SM observables whereas the light blue bands denote the solutions considering resonances in observables from Ref. [17] 
parametrization. The red error bars denote the solutions obtained using data (as highlighted in contours is Fig. 1). Including 
the resonances with the parametrization used in Ref. [17], the solutions for Pi and P 2 are unaltered and superimpose with the 
‘stars’ completely, (see text for details) 


affected by resonances. This is easily seen if we consider 
a situation where NP is absent and all the parameters for 
resonances (strength, phase etc.) are known, the observ¬ 
ables calculated using Eqs. (9)-(13) should agree with 
the experimental measured observables. Thus the con¬ 
sistent set of Eqs. (9)-(13) must provide the same set of 
parameters that we would have started with, as best fit 
solutions. In the absence of NP the measured observables 
should result in the solutions for parameters matching 
with SM values. Since Pi and P 2 are unaffected by reso¬ 
nances their best fit solutions also remain unaffected by 
it. Our best fit values of Pi and P 2 differ from the SM 
estimates and this discrepancy cannot be accounted for 
by resonances. 

To establish the above arguments we further un¬ 
dertake an extensive study illustrated in Fig. 3. We 
choose the region q 2 > 11 GeV 2 as resonance ef¬ 

fects can be dominant here and assume SM form 
factor values of the observables from lattice calcula¬ 
tions [14]. The solutions for Pi and P 2 are ob¬ 
tained using Eqs. (9)-(13) for seven different q 2 points; 
11 GeV 2 , 12 GeV 2 , 15 GeV 2 , 16 GeV 2 , 17 GeV 2 , 18 GeV 2 
and 19GeV 2 . The observables F_l, Fl , A 4 , H| B and 
A l are SM estimates calculated using lattice form fac¬ 
tors. These seven solutions of Pi and P 2 are denoted 
by ‘star’ symbols in the corresponding plots. The black 
central line with gray band is the form factor estimate 
(mean with ±1 a error) of Pi and P 2 . It can be seen 
that the set of Eqs. (9)-(13) are completely consistent 
with SM structure and produces expected solutions. In 
case the solutions were completely analytically obtained, 
the ‘stars’ should sit on the black curves. However the 
solutions for hadronic parameters are very complicated 
and has been evaluated numerically, resulting in small 
shifts that are visible. The blue error bars are the so¬ 


lutions for Pi and P 2 using the bin average values of 
SM observables. It can be seen that as the Eqs. (9)- 
(13) are valid at each q 2 point, bin averaging has in¬ 
duced some shifts in the solutions. However the results 
are in agreement within ±ler confidence level region. To 
illustrate the effect of resonances we have considered the 
parametrization from Ref. [17]. The charmonium bound 
states J/ijj(lS), 0(2S), 0(3770), -0(4040), 0(4160) and 
0(4415) are included in the mentioned five observables. 
Interestingly, the change in the value of observables in¬ 
cluding the resonances affected the solutions for £, u±_ 
and «||, however, solutions to Pi and P 2 remained unal¬ 
tered (upto second decimal place), hence, the solutions 
completely superimpose with the ‘stars’ obtained with¬ 
out resonance contributions. We have also investigated 
the effect of resonances in the bin average where the ob¬ 
servables are evaluated with lattice form factors including 
the above mentioned resonances and the solutions to Pi 
and P 2 are shown in light blue bands for the three q 2 
bins 11.0 < q 2 < 12.5 GeV 2 , 15 < q 2 < 17 GeV 2 and 
17 < g 2 < 19 GeV 2 . In this case the results with and 
without resonances do not completely superimpose how¬ 
ever are quite consistent within ±ler error bars. These 
solutions are also shown in Fig. 1 and 2, in same light 
blue bands, for the relevant bins where resonance effect 
may in principle be significant. The red error bars are 
the solutions for Pi and P 2 obtained from data (as dis¬ 
cussed and highlighted in contours is Fig. 1) that have 
been shown here again for convenience. We reiterate that 
effect of resonances in Pi, P 2 solutions is independent of 
the parametrization choice as the solutions do not depend 
on Wilson coefficient Cfj and our conclusions derived for 
Pi and P 2 parameters are unaffected by resonance effect. 
It is justified that bin average can induce some errors in 
the solutions. However, we have allowed a shift in the 
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observable values (in Eq. 23 and Table. IV) of more than 
the 1(7 error for each observable which hopefully is suffi¬ 
cient to compensate such effects. 

It is important to note that in our analysis no hadronic 
estimates are used to solve for the five parameters 
from exactly five measurements. Whereas, in other ap¬ 
proaches, when considering the same B —> K*££ mode all 
six form factors, Wilson coefficients and non factorisable 
corrections based on conservative estimations are needed. 
We compare Pi and P 2 obtained from experimental data 
alone, with the three form factor V, Ai and A 12 to which 
they are related as theoretical inputs. The form-factors 
Ti, T 2 and T 23 are not used in this comparison. Thus, our 
comparisons are different in nature and have reduced un¬ 
certainties, in terms of number of theoretical estimates. 
This may result in higher significance level of deviation 
observed here. 

The large q 2 region where the K* has low-recoil en¬ 
ergy has also been studied [18, 19] in a modified heavy 
quark effective theory framework which is a model inde¬ 
pendent approach. In this limit the number of indepen¬ 
dent hadronic form factors reduces to only three and one 
finds [3] that ro = r» = r± or equivalently uq = uy = u± 
must hold as long as non-factorizable charm loop contri¬ 
butions are negligible. We find that this relation does 
not hold for either of the bins 15 < q 2 < 17 GeV 2 or 
17 < q 2 < 19 GeV 2 . The values of uq, U\i and u± ob¬ 
tained from the fit with ±ler errors are listed in Table II. 
We note that u\ s receive problematic resonance contri¬ 
bution coming from Cg. To address this issue we have 
introduced more sytematics in measured observable than 
the one arising only from bin average effect. We have 
checked our analysis by doubling the systematics of the 
observables given in Table. IV of Appendix. B for the q 2 
range 11 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 and 
our results are stable with it. The actual significance of 
the deviations observed here can be obtained with the de¬ 
tailed study of resonance systematics which is a subject 
of an independent paper. However the significance level 
is evaluated by conservatively adding systematics vary¬ 
ing between 10% - 100% in the observables. The large 
discrepancies observed are equally hard to explain solely 
due to non-factorizable charm loop corrections and may 
be additional evidence of physics beyond the SM. 


q 2 range in GeV 2 

Mo 

u \\ 

Mi 

15 < q 2 < 17 

17 < q 2 < 19 

15 < q 2 < 19 

0.000 ±0.016 

0.166 ±0.014 

0.120 ±0.007 

0.013 ±0.153 

0.000 ±4.579 

0.004 ± 0.441 

0.367 ±0.025 

0.260 ± 0.048 

0.244 ±0.026 


TABLE II. The values of Mo, M|| and Mi obtained from fit 
to 311C 1 of LHCb data [4], In large q 2 region [18, 19] the 
equality Uo = it 11 = Uj_ is expected to hold if non-factorizable 
charm loop contributions are negligible. The errors in the 
value of My for the larger q 2 bin is unexpectedly large to draw 
any conclusions. Significant discrepancies which are too large 
to be solely due to non-factorizable charm loop corrections 
are observed between the values of u± and Mo in both bins. 


B. Testing relation between observables 


The relation between the observables for asymmetries 
A 4 , A 5 and Afb given in Eqs. (17) - (19) can also be 
tested using LHCb data [4]. In Fig. 4, top left panel, we 
compare theoretically calculated Aps mean values and 
± 1(7 errors (in yellow bands) with experimental mea¬ 
surements (red error bars) for the respective q 2 bins. 
All observables in the r.h.s of Eq. (19) (‘relation’ ) are 
assumed to be Gaussian distributions in data and the 
predictions for App in yellow bands are obtained using 
the expression of the ‘relation’. A very good agreement 
is evident for most q 2 regions, however, for the ranges 
11.0 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 a 
deviation of 2.1a and 1.8cr is observed. Similarly ‘re¬ 
lation’ for A 4 in Eq. (17) results in a very good agree¬ 
ment except for showing a discrepancy of 2.3a only in 
the 0.1 < q 2 < 0.98 GeV 2 bin, in right top panel of 
Fig. 4. The disagreement in the value of App and A 4 in 
some q 2 bins indicates that there is no set of form fac¬ 
tors and Wilson coefficients which can explain App and 
A 4 completely. Observables A 5 or equivalently Pg [20] 
are found to be in complete agreement i.e. within about 
± 1(7 deviation for all q 2 bins as shown in the two lower 
panels of Fig. 4. The solutions for A 5 and Afb have 
ambiguities. We chose the ambiguity for which the chi- 
squared deviations are the least. Our conclusions have 
no bearing on and do not rule out the observation made 
by LHCb in observable Pg in Refs. [4, 9]. The predictions 
of observable Pg derived from the relation is a signal of 
consistency of LHCb results. We note that the relation 
remains valid except in the presence of NP operators that 
result in modified new angular distribution. Hence we do 
not expect to see the discrepancy observed by LHCb [9] 
if right-handed currents or extra vector current such as 
Z' contributes to the decay. The discrepancy observed 
by LHCb depends on the comparison with model based 
calculation of form factors. Whereas, the predictions of 
these asymmetries made in this paper, are independent 
of any form factor values and depend purely on the gauge 
structure of SM. If the model dependent calculations of 
form factors are correct, signal of new physics may well 
be indicated in the bins suggested by Ref. [9]. We find 
that LHCb data indicates yet another independent dis¬ 
crepancy. 

IV. CONCLUSION 

In conclusion, we have used the 3 fb — 1 of LHCb data 
to determine some hadronic parameters governing the de¬ 
cay B —> K*£ + £~ assuming contributions from SM alone. 
We obtain the values of the form factors V(q 2 ), A\{q 2 ) 
and A 12 (g 2 ) that are used to describe the matrix element 
(iCIsy^PrAlP) directly from data. Very significant devi¬ 
ations are seen for the form factors V and Ai especially 
in the regions q 2 < 2.5 GeV 2 and q 2 > 6 GeV 2 . We 
point out that the presence of resonances in data can 
induce more systematic uncertainties in the fits. How- 










q 2 q 2 

FIG. 4. (color online) The mean values and ±ler uncertainty bands for asymmetries Afb, A 4 , Ag and Pg calculated using 
Eqs. (17) - (19) are shown in yellow, gray, green and brown bands, respectively. The error bars in red (dark) correspond 
to the LHCb measured [4] central values and errors for each observable for the respective q 2 bins. The predictions for the 
asymmetries are obtained using the relations among observables which are independent of any hadronic parameters and depend 
on experimental measurements of the other observables remaining in the corresponding relations. Sizable discrepancies are 
shown for Afb in 11.0 < q 2 < 12.5 GeV 2 and 15 < q 2 < 17 GeV 2 bins and for A 4 in the range 0.1 < q 2 < 0.98 GeV 2 . We 
note that the relations (Eqs. (17) - (19)) remain valid except in the presence of NP operators that result in modified angular 
distribution. Hence the presence of right-handed currents and any extra vector current such as Z' the relations will remain 
valid. 


ever in the view of absence of such a existence of res¬ 
onances in B —» K*£ + £~ data, we emphasize that the 
significant deviations observed in the form factor values 
can not be completely explained by resonances and non- 
factorizable contributions. We would like to point out 
that there exist major differences between the global fit 
approaches [ 6 ] to study the anomalies in b —> .s transi¬ 
tions and the approach adopted in our work. Our work 
relies only on B —> I\*££ decay mode, whereas, global fit 
techniques incorporate various decay modes and hence 
either use LCSR, Lattice based estimates of form fac¬ 
tors or treat form-factors as parameters in the fit pro¬ 
cedure. The number of inputs and fitted parameters 
differ making a number by number comparison of the 
different approaches difficult. Furthermore due to the 
absence of accurate estimates of non factorisable correc¬ 
tions, the global fit techniques rely on some conservative 
estimations of these corrections. However, the formalism 
we have developed parametrizes such corrections and the 
conclusions drawn here are independent of non factoris¬ 
able estimates. These are perhaps the reasons why we 
find larger significance. However, qualitatively we don’t 
see a significant disagreement with the other approaches 


as we do observe ~ 3<r discrepancy in Pi — P 2 plane in 
q 2 region [6 — 8 ] GeV 2 where observable P 5 also deviates 
by 2.7 a from its SM prediction. 


Further, a relation between form factors expected to 
hold in the large q 2 region as long as non-factorizable 
charm loop contributions are negligible, seems to fail. 
Finally, the relation between observables also indicates 
some deviations in the same regions where the form 
factors were found to disagree. The forward-backward 
asymmetry Afb deviates in the q 2 > 11 GeV 2 region, 
where as A 4 differs in the region q 2 < 0.98 GeV 2 . As 
the systematic error arises from the experimental mea¬ 
surements of observables in terms of binned dilepton in¬ 
variant mass are accounted, the magnitude of discrepan¬ 
cies observed would be hard to accommodate either as 
systematics from long distance resonance contributions 
or possible corrections to theoretical estimates. All these 
features can be understood if there are other unaccounted 
for operators contributing to the decay mode. In view of 
this, we speculate that these deviations are likely to be a 
signature of physics beyond SM. 
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Appendix A: Complex contribution ex estimates 
from data 


In Ref. [2] it was shown that the complex contributions 
£\ to the amplitude of the decay mode B -» K*i + t ~, can 
be taken into consideration. e\ can be solved in terms 
of iterative solutions proportional to the observables A 7 , 
A s , A 9 and a form factor ratio Pi. The expressions for 
all the three ex's are shown in Eqs. (76)-(78) of Ref. [2], 
They are reproduced here for convenience. 


£_L 


e l! 


Ed 


\/27rI/ 

A 9 P 1 A 8 P 2 A 7 PiP 2 r_i_ 

1 

0 

3 a/ 2 4 37rCio 

a/2 ?rl) 

Agr 0 i A 8 P 2 r|| A 7 P 2 r|| 

iS 

1 

0 

3\/2rj_ 4P rrx 37rCio 

x/2i rT/ 

AgPir 0 ! Agry A 7 Pir 0 

(r 0 -r ll )F± 

3V2P 2 r± 4r± 37rCio 


(Al) 

(A2) 

(A3) 


A point to be noted as explained in detail in Ref. [2], is 
that the (ea/I^/s are completely expressed in terms of 
observables and the form factor ratio Pi. However, these 
solutions are essentially iterative, since the r\ s and C 10 
are derived in terms of the primed observables that de¬ 
pend on e\. If (£\/V^ 2 ) are small as should be expected, 
accurate solutions for them can be found with a few iter¬ 
ations. In Ref. [2] the variation of £\ with P x was studied 
for 1 fb -1 LHC data and it was found that the solutions 
are not sensitive to the value of Pi. 

We generate a set of events for every bin, with each 
event consisting of randomly chosen values drawn from 
Gaussian distributions generated for each of the observ¬ 
ables Fl, Fj_, A 4 , A 5 . Afb, Ar, A s and Ag. The dis¬ 
tributions are generated using experimental results from 
Ref. [4], with the experimentally measured value as mean 
and the uncertainty as standard deviation. 


£\ are solved iteratively for every set of observables. 
We find converged iterative solutions for £\/y/Tf for each 
set of observables with the histograms shown in Fig. 5. 
The red (dark), light brown (lightest) and green his¬ 
tograms denote the solutions for £_l/-\/T}, £\\/\/l} and 
eg/y/Tf respectively for all the eight bins with q 2 range 
(0.1 - 0.98) GeV 2 , (1.1 - 2.5) GeV 2 , (2.5 - 4.0) GeV 2 , 



FIG. 5. (color online) The solutions for ex/\/l}, £\\/y/l } and 
£ 0 / x/Vf using distributions for first through eighth q 2 bins 
are depicted in red (dark), light brown (lightest) and green 
respectively. All the Ea/a/T/s are consistent with zero. 


(4 - 6) GeV 2 , (6 - 8) GeV 2 , (11.0 - 12.5) GeV 2 , (15 - 
17) GeV 2 and (17-19) GeV 2 . 

We have also quoted the mean and ±1 <t errors for 
each £\/x/Tf in Table. Ill calculated from the distribu¬ 
tions shown in Fig. 5. It can be easily seen that all the 
mean values of £a/\/I} are consistent with zero. From 
Eqs. (37)-(40) of Ref. [2], the contributions from imagi¬ 
nary part of the amplitude to the observables Fl, Fy, F± 
and A 4 are quadratic in the corresponding £\/x/^f and 
thus are negligible. 
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q 1 2 3 4 5 6 in GeV 2 

e±/Vff 

e ll/v^> 

£ o / v/l> 

0.1<(? 2 <0.98 

-0.048 ±0.116 

-0.047 ±0.103 

0.020 ±0.111 

1.1<q 2 <2.5 

-0.010 ±0.078 

-0.010 ±0.078 

0.078 ±0.172 

2.5<q 2 <4.0 

-0.009 ±0.079 

-0.008 ±0.080 

-0.025 ±0.212 

4-0<q 2 <6.0 

-0.026 ±0.097 

0.014 ±0.093 

0.032 ±0.234 

6.0<q 2 <8.0 

-0.011 ±0.088 

-0.046 ± 0.078 

-0.132 ±0.129 

11.0<q 2 <12.5 

-0.011 ±0.050 

0.038 ± 0.074 

-0.078 ±0.114 

15.0<g 2 <17.0 

-0.000 ±0.067 

-0.027 ±0.071 

0.020 ±0.072 

17.0<q 2 <19.0 

0.006 ± 0.076 

-0.090 ±0.090 

-0.040 ±0.088 


TABLE III. The E\/y/Tf mean values with ±1 <t errors from 
Fig. 5 


Appendix B: Systematic uncertainty evaluation 


q 2 range in GeV 2 

Fl 

F[ 

Al 

a 2S 

/ipB 

iOS 

Al 

0.1 < q 2 < 0.98 

0.014 

0.230 

0.088 

0.002 

0.016 

1.1 < q 2 < 2.5 

0.223 

0.151 

0.036 

0.034 

0.010 

2.5 < q 2 < 4.0 

0.164 

0.223 

0.064 

0.013 

0.004 

4.0 < q 2 < 6.0 

0.069 

0.138 

0.021 

0.002 

0.008 

6.0 < q 2 < 8.0 

0.132 

0.165 

0.028 

0.020 

0.019 

11.0 < q 2 < 12.5 

0.029 

0.063 

0.006 

0.051 

0.023 

15.0 < q 2 < 17.0 

0.019 

0.048 

0.027 

0.036 

0.023 

17.0 < q 2 < 19.0 

0.109 

0.020 

0.039 

0.077 

0.053 


TABLE IV. The systematic uncertainties for each observables 
F±, Fl, A 4 , Ap B and A 2 are shown. The values denote 
magnitude of maximum deviation of the bin average central 
value with the fitted q 2 polynomial within every q 2 bin. 

We discuss the evaluation of systematic uncertainties 


arising mainly due to bin average effect of observables. 
As written in Eq. (22), the shift O s in each observable is 
calculated for each q 2 bin, by considering the maximum 
deviation of the bin average value of the observable O 
from a fitted q 2 polynomial of entire range. It is high¬ 
lighted in Fig. 6 where red error bars are LHCb measure¬ 
ments and gray curves represent best fitted polynomial 
in q 2 for 14 bin LHCb data. We use 14 bin measurement 
(based on the method of moments [13]) from LHCb to 
fit the polynomial in q 2 , rather than the 8 bin data set 
as it provides more information to determine the shape 
of the polynomial for entire q 2 region. The black dashed 
line denotes the maximum deviation of bin average cen¬ 
tral value of the observable with the q 2 function for the 
region 6 < q 2 < 8 GeV 2 and O s is the length of the line 
for observable O. Similar technique is applied for other 
q 2 bins also and the values of systematic errors are given 
in Table. IV for all observables. 

It should be noted that as discussed in Sec. Ill finite 
lepton mass can affect the analysis in the first two q 2 
region namely q 2 < 2.5 GeV 2 and in the absence of a 
measurement of asymmetries A 10 and An [2] we have 
to rely on some hadronic estimates. This in principle 
may cause more uncertainties and we took a conservative 
approach by considering two times the O s values for all 
observables given in Table. IV for the two bins 0.1 < q 2 < 
0.98 GeV 2 and 1.1 < q 2 < 2.5 GeV 2 . 

We emphasize that resonances in our analysis will only 
affect the fitted function in q 2 , which in turn will induce 
more systematic uncertainties to the observables. We 
have checked the y 2 fit (in Sub-sec. Ill A) by increasing 
the systematic uncertainties two times of the values given 
in Table. IV for the regions 11 < q 2 < 12.5 GeV 2 and 
15 < q 2 < 17 GeV 2 and our results are stable with it. 
However a detailed study of resonance systematics on this 
decay mode is currently going on and will be a subject 
of an independent paper itself. 
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q 2 function for bin 6 < q 2 < 8 GeV 2 . The length of these black lines are denoted by Fj_, F[, A 4 , Ap B and Aj s , respectively. 
Similar lines can be drawn for other q 2 bins also and the values of systematic errors are given in Table. IV for all observables. 
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